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We examine the recent suggestion that V— and CP— odd effects in QCD matter can induce elec- 
tric charge asymmetry with respect to reaction plane in relativistic heavy ion collisions. General 
arguments are given which confirm that the angular momentum of QCD matter in the presence of 
non-zero topological charge should induce an electric field aligned along the axis of the angular mo- 
mentum. A simple formula relating the magnitude of charge asymmetry to the angular momentum 
and topological charge is derived. The expected asymmetry is amenable to experimental observation 
. . , at RHIC and LHC; we discuss the recent preliminary STAR result in light of our findings. Possible 

■ implications of charge separation phenomenon in cosmology and astrophysics are discussed as well. 
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I. INTRODUCTION 



I— !■ The strong CV problem remains one of the most outstanding puzzles of the Standard Model. The best known 
I solution of this problem is provided by the axion, see the original papers [H,[1|,[11j and recent reviews [1]. Thirty 
years after the axion was invented, it is still considered as a viable solution of the strong CV problem, and one of the 
O i' plausible dark matter candidates. At present, there are several experimental groups searching for the axions, see Q 
O and references therein. The axion scenario suggests that in the Early Universe V and CV invariances were strongly 
violated, while at present time they are preserved by strong interactions. V and CV conservation in strong interactions 
\ may thus be related to the vacuum structure. Moreover, recent progress in understanding the vacuum structure of 
CsJ ■ gauge theories, especially supersymmetric ones d, 0, Hi, points to the possible existence of V and CP-odd vacuum 
J> ' states. It was proposed that QCD vacuum in the vicinity of the deconfinement phase transition can also possess V— 
and C7^— odd metastable vacuum solutions [§] possibly inducing CP-odd correlations of the produced hadrons (loj 
in relativistic heavy ion collisions. Such "P-odd bubbles" are a particular realization of an excited vacuum domain 
which may be produced in heavy ion collisions [ll[ , and several other realizations have been proposed before (T^ . IT3| . 



\^ , The study of metastable V and CP-odd domains at RHIC provides a unique link to the QCD phase transition 
in the Early Universe when 9 parameter was likely not zero. Such a large CV violation during the QCD phase 

, transition might lead to a separation of matter and anti-matter in the Early Universe. As is well known, all terrestrial 

' experiments (conducted when 6 has already relaxed to zero, as the axion solution of the strong CV problem suggests) 

^ ' show very tiny CV violation which is not sufficient to explain the baryon number density observed today. We speculate 

• ^ in Section IV on a number of other astrophysical/cosmological implications which nonzero in a few moments after 

/\ the Big Bang might have. 

I Several dynamical scenarios for the decay of P-odd bubbles have been considered and numerical lattice 

calculations of the fluctuations of topological charge in classical Yang-Mills fields have been performed [H, ITg*! . A 
closely related topic is the role of topological effects in high-energy scattering and multi-particle production which 
has been studied in Refs . [l 8| . [ l9| . [20j . (2]| . The studies of P- and CP-odd correlations of pion momenta [13, HH, 
including those proposed in ref[24!|, have shown that such measurements are in principle feasible but would require 
large event samples. 



II. CHARGE SEPARATION EFFECT 



Recently, it was proposed that P and CP violation should manifest itself in heavy ion collisions through the 
separation of electric charge with respect to the reaction plane (255 . This phenomenon is a consequence of topological 
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charge fluctuations in the presence of non-zero angular momentum (and/or magnetic field) in off-central heavy ion 
collisions. The experimental detection of this effect is possible by the method proposed in ref [26j which correlates 
the emission angles of charged particles with the reaction plane on the event-by-event basis. First preliminary result 
of such a study has been presented recently by STAR Collaboration at RHIC [27|. The goal of this Letter is to 
re-examine the charge separation effect from a different point of view, to study the dynamics of this phenomenon in 
more detail and to provide numerical estimates needed for the interpretation of experimental results. 

Let us start by considering the effects of angular momentum on QCD matter. In experiment, hot QCD matter 
is produced by colliding relativistic heavy ion collisions, and such collisions occur at a finite impact parameter b. 
The initial angular momentum of the system is classically given by ^ and is huge at high energies ^/s; much 
of it is carried away by spectator nucleons but the produced hot matter nevertheless possesses significant angular 
momentum which combined with the gradients of pressure leads to the azimuthal anisotropy of hadron production. 
In the absence of parity violation, the angular distribution of the produced hadrons is symmetric with respect to 
the angular momentum axis (or equivalently, about the reaction plane) which of course fluctuates event-by-event. 
When expanded in harmonics, this angular distribution possesses contributions from even (as required by parity for a 
symmetric around zero rapidity interval) orbital momenta I. The harmonics of at least up to Z = 4 have been reliably 
identified in the data - see e.g. [1^ and references therein. 

What are the effects of the angular momentum on the constituents of matter? Larmor theorem states that classical 
angular precession with frequency fl for a particle of mass m and electric charge e is equivalent to the motion in 
magnetic field of strength 

B = (1) 

e 

Even more important for what follows is the fact that this correspondence persists on the quantum level; in particular, 
rotation of Dirac particle along the loop bounding an area S yields the phase shift of the wave function given by 



A<i>rot ^2m J n dE, (2) 
equivalent to the Aharonov-Bohm phase acquired by the motion in magnetic field B threading the same loop: 

A$B = e / B dS. (3) 



This phenomenon is known as Sagnac effect; it has been firmly established experimentally (see refs [29] for reviews). 
Having the correspondence ([1]) in mind, we will throughout the paper substitute the influence of angular momentum 
by the action of a background magnetic field. Such a substitution is natural also because off-central heavy ion collisions 
create magnetic fields aligned perpendicular to the reaction plane in the collision region occupied by hot matter. 

Let us assume, following the proposals outlined above, that a dynamical local fluctuation of 6'-angle can be excited 
in QCD matter. From the viewpoint of the effective lagrangian, 6{x) is equivalent to a pseudo-scalar flavor-singlet 
quark-anti-quark field which couples to electromagnetism through the electric charges of its quark constituents as 
prescribed by axial anomaly: 

'-¥'-IS'.N,^^^{±){e.b). (4) 

where the sum runs over quark flavors /, and N^, is the number of colors. Let us now minimize the action density Q 
with respect to the electric field E: 

^ = E + N,yi,-(^]B^O; (5) 
SE ^ 47r2 \NfJ ' ^ ' 

we see that the magnetic field in the presence of 6* ^ generates an electric field E ^ 9 ■ B. 

To understand this phenomenon better, we start by considering an example of a field B provided by a magnetic 
monopole of magnetic charge g: B = {g/r^) n. Using in ([5]) Nf = 1, A^c = 1, e/ = e as appropriate in this Abclian 
case, we get 



r2 V 27r y \2tt) ' 



n^-. (6) 

r 
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Replacing (e • g)/2-K = 1 (or in general an integer) in accord with Dirac's quantization condition in eq. (O we 
reproduce the famous result by Witten t30ll : magnetic monopole in the presence of 9 receives the electric charge equal 
to — (e 9)/2tt and becomes a "dyon" [311, [13], [l^l- In fact, dyons in QCD have been known to induce chirality 
violation and were proposed as a solution to the problem [3^ . [35| . A related phenomenon in grand unified 

theories is the proton decay induced by a magnetic monopole ~ Callan-Rubakov effect |3^ , [s^l • 

Let us now consider the case of a uniform magnetic field Bz pointing in the z direction; we will assume that the 
field does not depend on x and y coordinates, so we are dealing with an effectively 2-dimensional theory. Note that 
according to ([1]) this is the field which is equivalent to the action of angular momentum aligned along z. In this 
case, because of the quantization of the flux ([3]) required by the single-valucdncss of the particle wave function, we 
substitute J (Px^B^ = $/e = 2^1 /e into ([51) ; Hs an integer. We thus get (taking Nf = 1, = e ), 

o / e 



where L is the size of the system. Therefore, the electric field along z will be induced in the presence of nonzero 9^0 
when magnetic field is applied along the z direction. This is essentially the generalization of Witten's result [lO] for 
a non-monopole like magnetic field. The situation here resembles the 2d Schwinger model with 9^0 [s^ - indeed, 
this theory can be considered as occurring in a background electric field of strength (e 9)/2tt jssj . 

Now we want to use the analogy ([T]) to argue that an electric field will be induced even when an external magnetic 
field is zero B = 0, but instead the system is rotating with an angular velocity fl. In other words we anticipate a 
relation E ~ (%^) ^ which replaces eq. ([5]) when the magnetic field is replaced by V,. Amazingly, such an exact 
relation indeed can be derived using the anomalous effective lagraiigian approach [39| in the presence of matter with 
nonzero chemical potential fj,, and its generalization to nonzero 9 [40]. 

One can show that the electric field will be induced for the rotating system according to the following 

relation, see Appendix for details. 

where /x/ is the chemical potential of the species /. This result is quite remarkable - it indicates that the angular 
momentum in the presence of ^ generates a background electric field even in the absence of magnetic field. 

The electric field ([8]) aligned along the vector of angular momentum (and thus perpendicular to the reaction plane) 
will act on charged particles and cause the electric charge separation signaling the violation of V and CV invariances 
arising from 6* 7^ 0. One can now easily estimate the magnitude of the expected effect once the crucial ingredient (|8]) 
is known. A detailed modeling of the charge separation phenomenon would require the knowledge of charge transport 
properties, the geometry of the system, etc. Fortunately, we do not have to do all these additional assumptions 
because an exact formula (for charge density accumulated at the opposite sides z = ±L/2 of the system of size L) 
can be derived, see Appendix, 



where we assume that the system containing 6* ^ domains (7^— odd bubbles) has sufficient time to equilibrate^. 

This result again indicates that the angular momentum in the presence of ^ 7^ generates a background electric 
field even in the absence of magnetic field. Moreover, the magnitude of the charge density does not depend on the size 
of the of the system along z direction. Of course, this result is a consequence of topological (anomalous) nature of the 
phenomenon. Now the result Q is very easy to interpret: the induced charge separation ([9]) will result in induced 
z independent electric field ([8]) as long as the system can be considered effectively two-dimensional, Lx^Ly ^ Lz. 
As is well known, the electric field Ez between two infinitely large charged plates with charge density axy is equal 
Ez — cTxy Our equations ([8|) and ([9]) satisfy this equation for an infinitely large capacitor. 



^ Some arguments indeed do suggest that the system has sufficient time for the 8 vacua to be formed at RHIC [l4j . 
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FIG. 1: Preliminary STAR data from 27j. Charged particle asymmetry parameters as a a function of centrality bins selected 
on the basis of charged particle multiplicity in |j7| < 0.5 region. Points are STAR preliminary data for Au+Au at -^sjvjv = 62 
GeV: circles are a^., triangles are and squares are a+a_, see text for details. Black lines are theoretical prediction [5^ 
corresponding to the topological charge \Q\ = 1. 



III. NUMERICAL ESTIMATES 



We now want to get from ([9]) an estimation for an excess of quarks over antiquarks, say, in the upper hemisphere, 
i.e. above the reaction plane (the lower hemisphere will obviously have an equal excess of antiquarks over quarks). 
Therefore, integrating ([S]), see eq. (|22p in Appendix, and assuming /i„ — /i^ we have 



(10) 



For a numerical estimate we put Nc — 3, Nf = 2 to arrive at 



Nq^q = 2g ~ 3e 



(11) 



which shows that the asymmetry vanishes for strictly central collisions when the angular momentum is equal to zero, 
I — and increases for more peripheral collisions. Note also that this formula and Eq.([7|) clearly resemble the 2d 
Schwinger model with 9^0, see [33|. 

In the presence of a collective expansion of the produced system (as evidenced by the present RHIC data), the 
spatial separation of charges will result in the momentum anisotropy of the emission of positive and negative particles. 
Depending on the sign of 9 and I which fluctuate on the event-by-event basis, an upper hemisphere can thus have 
either excess of quarks over anti-quarks or viceversa. Taking for the sake of estimate ~ 1, using ^ = 4 for a 
semi-central Au — Au event as identified at RHIC ^281] , and e = \/47ra^, we find that a typical event at RHIC would 
have on the order of 5 more quarks than antiquarks in the upper hemisphere (above the reaction plane) , and an equal 
excess of antiquarks over quarks in the lower hemisphere. Assuming, in accord with the local parton-hadron duality, 
that the multiplicity of quarks and antiquarks at hadronization is approximately equal to the multiplicity of produced 
hadrons^ which for semi-central collisions at RHIC is 0(500), wc find an asymmetry between quarks and antiquarks 



A = 



N, 



0(3-5%), 



(12) 



q+q 



which agrees with the estimate of ref "25"] . The effect of this magnitude may be amenable to experimental observation; 
a precise way for extracting the charge asymmetry is offered by the mixed harmonics method as proposed by Voloshin 



^ It is assumed that the produced gluons split into quark-antiquark pairs before they hadronize. 
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[201 . In this approach, one studies the correlation of the produced particles with the reaction plane on the event-by- 
event basis. The first preliminary studies of charge separation effect have been performed by the STAR Collaboration 
[27I I at RHIC^. The parameters a± plotted in Fig.l in the absence of directed flow (which should be the case for a 
rapidity interval symmetric around rapidity y = in a symmetric collision) can be represented as 



\ k.l 



^PRp)sm{d-ijRp)), (13) 



where V'-RP is the reaction plane. In the absence of dynamical T'-odd correlations in the emission of positive and 
negative particles the average factorizes into the product of the averages of two sines of the emission angle, which 
should vanish because the system is symmetric w.r.t. the reflection around the reaction plane. The (unbroken) charge 
symmetry implies a+a+ = a_a_ > 0; because the negative and positive charges in the presence of charge separation 
effect move to the opposite sides of the reaction plane, we also expect a+a- < 0. Assuming the angular distribution 
of the produced particles to be [2^ 

dN 

— - l + 2asin(/), (14) 
dq) 

where (j) is the emission angle w.r.t. the reaction plane, the asymmetry A which we have estimated is related to 
the parameter a by a simple relation: A = 7ra/4. The magnitude of the correlation presented in Fig.l thus appears 
comparable to the expected one. The vanishing of the asymmetry in central collisions (with angular momentum 
I = 0) and its increase towards more peripheral ones is also qualitatively consistent with expectations. We therefore 
urge the continuation of experimental studies; we also note that a more direct way of observing the quark-antiquark 
asymmetry would be provided by the studies of baryon-antibaryon production. 

We would like to complete the section with a cautionary remark on the accuracy of our estimates: all theoretical 
formulae discussed above have been derived at zero temperature using the effective lagrangian approach. We expect 
that both the absolute value of the asymmetry and the real time dynamics of charge separation will change in a 
more realistic calculation at T ^ 0. However, since the effect we consider is of topological nature, it should persist at 
T 7^ 0. In addition, the signal can be reduced if a large number of 6 domains is produced in each event; however we 
do not think this is likely because the typical size of the domain is not much smaller than the size of the produced 
matter. A more detailed treatment will be presented elsewhere; the main goal of the present letter is to demonstrate 
the topological nature of the charge separation effect as a manifestation of strong CV violation in QCD. 



IV. POSSIBLE IMPLICATIONS FOR ASTROPHYSICS AND COSMOLOGY 



A. Little Bang versus Big Bang 



• In the concluding section of this paper we would like to speculate on the possible consequences of the charge 
separation phenomenon for the evolution of the Universe shortly after the Big Bang during the QCD phase transition 
at To ^ 170 MeV. In the axion framework, see refs [l| - it is generally assumed that parameter 9 was not zero when 
the temperature of the Universe T was larger than Tq, T > Tq. The 6 angle starts to relax to zero only at the time of 
the QCD phase transition as a result of formation of the chiral condensate: the axion acquires its mass ^ mq{qq) 
while the axion potential Va{9) being almost flat at T > Tg develops the global minimum at 6* = when T slowly 
decreases to T < Tq due to the expansion of the Universe. Eventually, the axion field relaxes to the minimum of the 
potential at = such that strong CV problem is resolved at present time no matter what the original magnitude of 
6 before the QCD phase transition was. This picture represents the standard resolution of strong CV problem, see the 
original papers 2],[3], and recent reviews 4]. However, we stress that while the parameter 6 vanishes at present 
epoch, = 0, it was of order one, ~ 1, during the QCD phase transition. We accept this standard picture of the 
axion evolution for the discussions which follow. In different words, we treat 9 as being a nonzero parameter of order 
one during the QCD phase transition time when the charge separation phenomena (the subject of the present work) 
may occur. The crucial ingredient of this picture is the existence of strong CV violation during that time. In particular 
this picture does imply that quarks and anti-quarks interacted very differently at that time in a stark contrast with 



^ We stress the preliminary nature of the reported data and the articulated by the authors of ref [27| necessity to analyze all possible 
systematic errors. 
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the results of present day experiments which show that CV symmetry is preserved with an extraordinary accuracy. 
This new source of very large CV violation in the Universe at time when T ^ Tq may be a crucial missing element for 
solving the puzzle of baryogenesis; we will comment on this below. 

• Another important ingredient which is always present when the physics of 6 angle is discussed is the existence of 
domain walls which form during the QCD phase transition, see Sikivie in [4] for a review on the subject. The existence 
of these domain walls is a direct consequence of discrete symmetry 9—^6 + 2?™ related to the angular nature of the 
9 parameter. The axion domain wall network forms at the QCD phase transition and eventually decays by emitting 
axions as well as by forming the closed surfaces, bubbles which will also eventually decay to axions if nothing stops 
this violent collapse. However, if the number of quarks trapped inside of a bubble is sufficiently large, the collapse 
stops due to the internal Fermi pressure, see [4l[ for detail discussions. Compact objects containing large numbers 
of (anti)quarks can thus be formed, which could play a role of cold dark matter candidates. The separation of 
charges which may happen during this period of evolution of the Universe is absolutely necessary for this scenario 
to work. Indeed, in this case the excess of quarks over anti quarks (or vice versa) inside the bubble with definite 
baryon charge could be macroscopically large preventing the bubble from collapse. The 9 values on different sides of 
the axion domain wall separating internal and external regions of the bubble are different. This feature models the 
domain which may be formed at RHIC with 9^0 inside the domain and 9 — outside of it. An important difference 
however is that in the Early Universe the domain size is determined by the (inverse) axion mass and is large, 
whereas at RHIC a much smaller domain size is expected. 

• One more crucial element of the present work is demonstration that CV odd combination {9E -Q) 7^ is correlated 
on the scale of the fireball size, large compared to AgJ^^. Similar CV odd correlation {9E ■ B) ^ is expected to 

occur in the early Universe where fl in our calculations (relevant for RHIC physics) is replaced by magnetic field B 
presumably induced on macroscopical scales in the early Universe'*, see [ol. [42l|. 



B. Separation of charges: consequences for cosmology 

•The objects with a large number of trapped (anti)quarks, if formed, would be absolutely stable macroscopical 
objects with the size of order m~*. They may serve as the cold dark matter candidates [4l|. The nuggets despite 
being baryonic configurations would behave like nonbaryonic dark matter objects because of tiny interaction (per unit 
mass) with the visible baryons, see [i^ for details. The stability of the nugget is due to the fact that the mass of 
the nugget is smaller than the mass of a collection of free sep arated nucleons with the same baryon charge. In many 



aspects these nuggets are similar to Witten's strangelets [4j| with very similar phenomenology. 



If the separation of matter from antimatter did indeed occur on macroscopical scales during the QCD phase 
transition, it may lead to a number of profound consequences. In particular, we speculate that two of the outstanding 
cosmological mysteries - the nature of dark matter and baryogenesis - may be explained by the idea that dark matter 
consists of dense nuggets of matter as well as antimatter. Due to the separation of matter from antimatter, these 
nuggets may form at the same QCD phase transition where "conventional" baryons (neutrons and protons) form. 
This would provide a natural explanation for the similarity of the dark matter and baryon contributions to the critical 
mass, floM ~ Sr^B, and for the baryon asymmetry in our part of the Universe through a charge separation mechanism 
driven by the strong CV violation due to the 9 term in QCD (see [1^ for details). Indeed, an "excess" antimatter in 
this case is locked away in antimatter nuggets; no additional baryon number violation is thus required to explain the 
observed matter/antimatter asymmetry in our part of the Universe. 

It is important to remark here that for this mechanism to be operative there is no need for any fine tunings (which 
arc an unavoidable ingredient of most baryogenesis scenarios); any efficiency of the nugget formation larger than 10^^*^ 
would be sufficient to explain a similarity of two contributions, ^dm ^ ^b- This is due to the fact that the present 
day baryon density/entropy ratio ub/s ^ 10^^° is 10 orders of magnitude smaller than it was at the QCD phase 
transition epoch due to the quark-antiquark annihilations. Therefore any quark-antiquark asymmetry (manifesting 
itself through the difference in the number and/or size of nuggets and anti- nuggets) will survive until the freeze-out 
time when quark-antiquark annihilation stops. A large CV violation due to the 9 term at early times provides a 
difference of order one between number density of nuggets and anti-nuggets such that a symmetric picture with an 
equal number of nuggets and anti nuggets cannot be realized^. Baryon charge conservation will automatically provide 



The scales relevant for the axion physics where correlation {6E • S) ^ is expected to occur are the scales at which the 8 parameter 
varies, i.e. on the scales of order of where ttIq is the axion mass, ~ {10~^ — 10~^) eV. 

^ Again, there is no fine tuning requirements. The only requirement is the following: CV violation must be sufficiently large > 10^^" at 
the time of the QCD phase transition. 
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an appropriate baryon number density to be of order of Qdm irrespectively of any initial conditions or specific 
details of the nugget formation mechanism as long as it is sufficiently effective, i.e. occurs with a relative probability 
> 10^^". Of course, all these attractive features of our scenario are related to the fact that both contributions 
originate at the same time at the same QCD phase transition and are made of the same constituents, quarks, see 
the original paper for the details. One should note that in general the relation ~ ^dm, within one order of 
magnitude, between the two different contributions to 17 is extremely difficult to explain in models that invoke dark 
matter candidates not related to the ordinary quark/baryon degrees of freedom. 

• Such a "counterintuitive" proposal does not in fact contradict an y of the known observational constraints on dark 
matter or antimatter in our Universe due to the three main reasons 1431: 



1) The nuggets carry a huge (anti)baryon charge \B\ « 10^° - 10^'^, so they have a macroscopic size (m^^ ^ 10^'^cm) 
and a tiny number density; 

2) They have nuclear densities in the bulk, so their interaction cross-section per unit mass is small a/M « 10"^'^ - 
10~^ cm^/g. This small factor effectively replaces a condition on weakness of interaction of conventional dark matter 
candidates such as weakly interacting massive particles (WIMPs); 

3) They have a large binding energy (gap A « 100 MeV) such that baryons in the nuggets are not available to 
participate in the Big Bang nucleosynthesis (BBN) at T sa 1 MeV. 

• In fact, this picture may also naturally explain a number of other mysteries in astrophysics and cosmology as we 
will now discuss. Several independent observations of the Galactic core suggest hitherto unexplained sources of energy 
and have detected an excess flux of photons across a broad range of energies. If charge separation indeed took place 
during the QCD phase transition, we may be seeing the consequences of this at present by observing the results of 
rare events of annihilation between the visible matter (electrons and protons) and antimatter nuggets. In particular, 
the observations which are very difficult to explain by conventional means include: 

1) SPI/INTEGRAL observations of the galactic centre have detected an excess of 511 keV gamma rays resulting 
from low momentum electron-positron annihilations. The observed intensity is a mystery: after accounting for all 
known positron sources, only a small fraction of the emission may be explained [46]. 

2) Detection by the CHANDRA satellite of diffuse X-ray emission from across the galactic bulge provides a puzzling 
picture: after subtracting the known X-ray sources one finds a residual diffuse thermal X-ray emission consistent with 
a two-tempcrature plasma with the hot component temperature close to T ~ 8 keV. According to the analysis [i^l 
the hot component is very difhcult to understand within the standard picture. Such a plasma would be too hot to be 
bound to the galactic center. The authors of ref. [47] also remark that the energy required to sustain a plasma of this 
temperature corresponds to the entire kinetic energy of one supernova every 3000 yr, which is unreasonably high. 

3) The flux of gamma rays in the 1-20 MeV range measured by COMPTEL represents yet another mystery. As 
discussed in [i^ the best fit models for diffuse galactic 7 rays fit the observed spectrum well for a very broad range 
of energies, 20 MeV- 100 GeV. It also gives a good representation of the latitude distribution of the emission from 
the plane to the poles, and of the longitudinal distribution. However, the model fails to explain the excess in the 1-20 
MeV range observed by COMPTEL in the inner part of the galaxy. So, what is the origin of the sources contributing 
to the energy band, 1 MeV < fc < 20 MeV ? 

It is quite remarkable that all these mysteries in different energy bands can be naturally explained if charge 
separation indeed took place in the Early Universe: an excess of the emission we observe now would be a result of 
annihilation between the visible matter and the nuggets of antimatter. Dark antimatter nuggets would provide an 
unlimited source of positrons (e^) within this framework as suggested in [49]. The resonance formation of positronium 
between impinging galactic electrons (e~) and positrons (e+) from the DM nuggets, and their subsequent decay, lead 
to the 511 keV line. Non-resonance direct e+e" 2j annihilation would produce [50| a broad spectrum at 1 MeV 
< fc < 20 MeV which is identified with the excess observed by COMPTEL. It is quite remarkable that another (also 
naively unrelated) puzzle, the diffuse X-ray emission observed by CHANDRA may also have a common origin with 
the 511 keV line and excess MeV radiation as argued in [5lj. This emission is due to proton (rather than electron) 
annihilation events releasing about 2mp ~ 2 GeV of energy per event. It occurs inside the nuggets, therefore, the 
available energy 2 GeV energy will not be radiated in the form of GeV photons, rather it will be radiated as X-rays 
in keV band [5l|. Other implications, in particular, may include generating the seeds of a primordial magnetic field 
at the QCD phase transition scale [ol. |4^. Is^ . 

• To conclude: several independent observations of the Galactic core suggest unexplained sources of energy. We 
speculate that the excess of photons may come from rare events of annihilation of visible matter with " dark" antimatter 
nuggets. The charge separation phenomenon which is the subject of this paper may induce the separation of matter 
and antimatter and the formation of nuggets during the QCD phase transition if 6* ^ at that time. Charge 
separation effect may thus hold the key to many mysteries of the Universe, and we urge the experimental studies of 
this phenomenon in heavy ion collisions at RHIC and LHC. 
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Appendix 

The main goal of this Appendix it to show how the eqs. ([l]),® follow from the anomalous effective lagrangian in 
the presence of nonzero chemical potential and nonzero 9{x). The corresponding physics of the anomalous effective 
lagrangian in the presence of 9 and fi is quite rich and deep with a number of connections to different branches of 
physics. In the present paper we shall give only a short sketch of how such relations arise without much details. We 
assume that 9{x) depends on coordinates and time even though for the purposes of the present paper we do not need 
to assume the existence of the axion. The main reason for that assumption is that we expect that a small domain 
with 6* 7^ can be formed at RHIC. This region will be treated as a macroscopically large region inside of a much 
larger region (representing our world) with 9 = 0. The 9{x) will be treated as a classical background external field 
without dynamical kinetic term. Nontrivial boundary conditions between the regions of 9 = and 9 ^ Q will play the 
crucial role in our discussions. The corresponding technique with nonzero ^ was developed in [39| and generalized for 
the axion field (nonzero 9{x)) in [ioj. 

We limit ourselves only to one term in the anomalous effective lagrangian which plays the crucial role for the present 
study, see [39] and [40] for details: 

Lg,v = -NcJ2i^ ■ ''"'^^WxV,)A,, (15) 
f ^ 

where is a fictitious vector gauge field introduced in (39j . Then the coupling of quarks to and to the usual 
electromagnetic gauge field are almost identical; the only difference is in the coupling constants. This new term 
(jlSp leads to a few very unusual phenomena. First, it leads to the separation of charge. Indeed, by definition, the 
local charge density is defined as 

^'_L^^^^^^. ,,.smv.) ^ A-.E ^ ^ (V" ^ 0) . (16) 

where 0, is defined as 2eijki^k = (diVj — djVi) and is nothing but the angular velocity of the rotating system. The result 
(|16p is very suggestive and implies that there is a charge separation phenomenon in the presence of 9 in the rotating 
system along the vector of the rotation fi. Indeed, assuming that f2 is directed along z we arrive to the following 
expression for the density charge Q/Hxy ( per surface area T,xy in the plane perpendicular to fi^) accumulated at 
z = ±L/2 where there is a transition between induced 9 state and vacuum^, 

/O rL/2 



l^xy Jo J 27r^ A// 

where is a value of 9{x ~ 0) inside the bubble for a given event while 9 ~ in the vacuum, outside the region of 
interest. A charge with the opposite sign will be induced at z = ~L/2. This is precisely eq. ([SJ used in the text. 
Let us rewrite eq. (fT5|) in the following way, 

L = lE^ + N,y^^-9(E-n), (18) 
f ^ 

where we added the standard kinetic term ■^E^. Without the anomalous term the ground state corresponds to 
(E) = 0. However, in the presence of the 9 term the minimum of the free energy corresponds to the state where 



here and in what follows we assume a cylindrical geometry defined by the length L along z and area Si 
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{E) 7^ 0. Indeed, minimization of (fT8|) with respect to E gives, 

which is precisely eq. ([8]) used in the text. As is known, the electric field i?^ between two infinitely large charged 
plates with charge density axy is exactly equal Ez — <Jxy- Our equations ()17|) and (I19p satisfy this famous equation 
for an infinitely large capacitor even though they were derived in a very different way. 

One more remark on eq. (|19p . An additional term in the lagrangian (jlSp is a total derivative, and therefore does 
not change the equations of motion. However, in order to quantize the system ()18p in canonical way, we have to 
express everything in terms of the generalized momentum tt^ — jj- which is defined as 




In the ground state we require that {iTi) — which is precisely the condition p9|l derived earlier. In different 
words, while the 9 term is a total derivative and does not change the equations of motion, it does influence the physics 
(separation of charges, induced electric field) in the presence of topological background field, flk ^ e.ijk{diVj—djVi) ^ 0. 

In order to integrate over the surface (Hlxy we have to understand the way how the rotation is quantized. The 
answer lies in the definition of the fictitious field = which is the local 4— velocity of matter with chemical 
potential /i/, see [s^jiO]- Then the coupling of quarks to and to the usual electromagnetic gauge field A^j_ takes 
the form: 

^ = - e/A,)^/7.V'/ (21) 

/ 

With this definition magnetic flux quantization as well as rotation quantization for a single species are very similar 
and take the form, 

j eB- (E^ 2-111 (22) 
y" Atf2 • (E = 2itI (23) 
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